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The multipole moments method is not only an aid to understand the deformation 
of the space-time, but also an effective tool to solve the approximate solutions of 
the Einstein field equation. However, The usual multipole moments are recursively 
defined by a sequence of symmetric and trace-free tensors, which are inconvenient 
for practical resolution. In this paper, we develop a simple procedure to generate the 
series solutions, and propose a method to identify the free parameters by taking the 
Schwarzschild metric as a standard ruler. Some well known examples are analyzed 
and compared with the series solutions. 
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I. INTRUDUCTION 

The asymptotically flat vacuum solutions with axisymmetry play an important role in 
general relativity. To understand the global structure of the metric caused by stars, we have 
the exact solution families of Schwarzschild, Curzon and Kerr metrics and some of their 
extensions to the electrovacuum solutions such as Reissner-Nordstrom and Kerr-Newman 
metrics [l]-[4]. After more than ten years of work, Manko et al. got an axisymmetric 
solution with five free parameters [5]. This solution is suitable for a rapidly rotating neutron 
star, in which the strong magnetic field and the oblate shape should be taken into account [6]. 

When getting a realistic exact solution becomes more and more difficult, the effective ap- 
proximate method is a good alternation. The multipole moments expansion was introduced 
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to interpret the physical meaning of the solutions to the Einstein equation. Geroch defined 
a recursive relation of multipole moments for the static asymptotically flat vacuum with 
axisymmetry by means of the timelike Killing vector [7, 8]. The norm of the Killing vector 
satisfies the Laplace equation in the 3-dimensional hypersurface orthogonal to the Killing 
vector, which is reduced from the Einstein equation. For the asymptotically flat space-time, 
the norm of the Killing vector and the metric is analytic near the spatial infinity, so the 
solution can be expanded in series of r -1 with multipole moments. The multipole moments 
can be recursively defined by a sequence of symmetric and trace-free tensors. This idea was 
generalized to the stationary case by Hansen in [9], where two different sets of multipole 
moments are employed. One is the mass moments and the other is the angular momentum 
moments. The two sets of moments can be generated from two potential functions defined 
by the norm and the twist of the time-like Killing vector. 

The formalism proposed in [8, 9] is in the covariant form. Some theorems are proved 
based on this formalism. In [10], Xanthopoulos proved that a stationary space-time is static 
if and only if all of its current moments vanish, and static space-time is flat if and only if all 
of its mass moments vanish. Beig-Simon [11] and Kundu[12] proved that, two spacetimes 
with the same multipole moments have the same space-time geometry at large radii, where 
the multipole expansion of the metric converges. This formalism was realized in the Weyl- 
Lewis-Papapetrou canonical metric[13, 14, 15]. 

An another formalism of multipole moments for realistic calculation of slowly changing 
systems or precisely stationary ones was reviewed and developed by Thorne in [16]. This 
formalism is more manifest for practical resolution, which is associated with a special kind of 
coordinate system called 'Asymptotical Cartesian and Mass Centered' (ACMC) coordinate 
system. In such ACMC-iV coordinate system, the metric can be expressed by multipole 
moments expansion, and the coefficients of r~(' +1 )(7 < N) terms are linearly combinations 
of the spherical harmonics of order I, I — 1, • ■ • ,0. The de Donder coordinate system belongs 
to ACMC-oo. The comparison and the relation between the two formalisms were discussed 
in [17]. This formalism in de Donder coordinate system and quotient harmonic one was 
developed into an effective tool for solving the approximate global solutions to the Einstein 
equation[18]-[20]. 

However, most of the previous works concerned mainly theoretical aspect of the multipole 
moments, which are very inconvenient for practical resolution of the Einstein's field equation. 
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Besides, the interpretation of the physical meanings of the parameters is still a general 
problem in general relativity [21, 22], although the primary motivation of Geroch to introduce 
the multipole moments was to clarify this problem. 

In this paper, we consider how to effectively solve the series solutions to the stationary, 
asymptotically flat vacuum with axisymmetry. This case is of important physical significant. 
At first, we show that the canonical form of Weyl-Lewis-Papapetrou coordinate system is 
an ACMC-oo one. Then we solve the series solution straightforwardly in this coordinate 
system. To understand the meanings of the free parameters in the solution, we set the 
Schwarzschild solution as a standard ruler. By comparison with this solution, we find that 
the bigger the absolute value of the dimensionless free coefficients, the larger the deformation 
and convection of the gravitating source. Noticing its generality, the series solution is an 
effective alternation to the exact one, and is helpful to understand the abundant structure 
of the space-time. 

II. THE SERIES SOLUTION IN WEYL-LEWIS-PAPAPETROU METRIC 

To describe the axisymmetric space-time, the canonical form of the Weyl-Lewis- 
Papapetrou metric is the simplest one[l]. The line element is equivalent to 

ds 2 = U(dt + Wdcp) 2 - V(dp 2 + dz 2 ) - U~ l p 2 d<p 2 , (2.1) 

where (t, p, z, tp) is the coordinates of the geometrical meaning near the cylindrical coordi- 
nates in Minkowski space. (U, V, W) only depend on (p, z). In this paper we take G = c = 1 
as units. Since the following calculations take the dimensionless form, the units will not be 
confused. The canonical form (2.1) has an important property, that is, 

Lemma 1. Except for the translation t = t + to,ip = ip + ip related to two Killing vectors 
{dtjdp) and z = ±(z + zq), the form of metric (2.1) removes other uncertainty caused by 
coordinate condition. 

This property can be checked as follows. The transformation (p, z) — > (p, T) keeping 
V(dp 2 + dz 2 ) = V{d'p 2 + dz 2 ) is a conformal transformation 

p = X[f(p±zi)], z = $[f(p±Zi)}, (2.2) 

where / is any given analytic function. But the restrictions p 2 dtp 2 = fPdip 2 and p > give 
/ = p ± (z + Zq)z, then we have z = ±(z + z ). 
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Except for the axial symmetry, the celestial body in equilibrium and the related stationary 
metric still have another top-bottom reflection symmetry, which means the metric will be 
the even functions of z by setting the origin (p = 0, z — 0) at the mass center. In what 
follows we only discuss this case, that is, (U, V, W) are even functions of z, although the 
discussion is also valid in the case including odd terms of z. So if we set the direction of z, 
the relation between coordinate and metric is completely fixed. 

The form (2.1) is not convenient for the following calculation, so we make a polar coor- 
dinate transformation 

p = rsin#, z = rcos9, (2.3) 

then the line element (2.1) becomes 

ds 2 = U{dt + Wd V ) 2 - V{dr 2 + r 2 d0 2 ) - U~ l r 2 sin 2 6d(p 2 . (2.4) 

Calculating the Ricci tensor R^ u , we get the independent components among the Einstein 
equations for the axisymmetric vacuum[l] 

St, + \ v + ^ dtU + _ r™!I + q^El s _ 2VRtt = o, (2 . 5) 

r r z r A U r z sin 9 

d 2 W - + ^d 2 e W + |(W ■ W) = ^(WR tt + R tip ) = 0, (2.6) 

where V = (d r , r _1 <%). And then V can be integrated from 

d r V d r U rsm 2 6 rrt2 2cot0 Qrr „ rr 1 rA2 \ 

ir + 17 {(W + —^drUdeU - - 2 (d e u) 2 j + 

((d r W) 2 + ^-drWdeW - - 2 (d e W 
It \ r r 

— sin 2 9(r 2 R rr - R ee ) - 2 cos 9 sin 9R r8 = 0, (2.7) 



r 



V U 2U 2 V r r 

^-(cot9[(d r W) 2 -\(d e W) 2 ]--d r Wd e w) = 

cos 9 sin 9{r 2 R rr - R ee ) - 2r sin 2 9R r0 = 0. (2.8) 

Substituting R^ = niT^ — ^g^T®) into the above equations, we get the dynamical equa- 
tions for interior metric, which can also be solved by series. 



5 



The solutions to (2.5)-(2.8) of most physical interests are the asymptotically flat cases, 
that is, the solution satisfies the boundary conditions at r — > oo 

U ^i_^ + 0(r" 2 ), W ^— sin 2 9 + 0(r~ 2 ), V^l + — + 0(r~ 2 ), (2.9) 

ry iy 

in which m is the total gravitational mass, and J the angular momentum of the source. 
(2.9) means r = oo is the regular point of the solution. On the other hand, the main part 
of derivatives in (2.5) and (2.6) are the linear elliptic operators. So the vacuum solutions 
(U,W) are analytic functions in the neighborhood of r = oo[ll, 12], and can be expressed 
as Taylor series of r , then we generally have 

^ r r 2 ^ ^ r n ' ^ ^ r n ' (2.10) 



u 



n=0 n=l 



where (R s , A n , B n ) are functions of ^ to be determined. For the regular metric, (U,V) are 
bounded even functions of z, or equivalently, (R s , A n , B n ) are bounded functions of cos 2 9. 

Substituting (2.10) into (2.5), and then expand it into Taylor series of r _1 , we can get 
the equations for the coefficients 

R" s + cot 9R' S = 0, (2.11) 
A% + cot 9A' + 2A = (R' s f + R 2 s . (2.12) 

The bounded solution of reflecting symmetry with respect to z = or 9 = |, we call it 
regular solution in what follows, is given by 

R s = 2m, A Q = l -R 2 sl (2.13) 

where m is a constant, which means the total mass of the source, m = corresponds to 
the Papapetrou class of solutions or flat space-time. The case of m < is unphysical. So 
only the case with positive mass m > has the most physical interest, and we only discuss 
this case in what follows. Obviously, R s is the corresponding Schwarzschild radius with 
dimension of length. 

By virtue of the excellent structure of (2.5) and (2.6), we have the following results. 
Theorem 2. For the regular solution (2.10) with starting value (2.13), it takes the 
following form 

Re 1 / Rq \ ( Re \ T ( R<i 



°° ~ ( R\ n 

W = R s J2 B n[ — ) , (2-15) 

n=l V r / 
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in which 

A n = a nk cos(2k6) , B n = b nk [cos(2k9) - 1] , (2.16) 

k=0 k=l 

where (d n k,b n k) are dimensionless constants to be determined, and N n = [^rp-] stands for 
the integer part of |(n + 1), namely, N 2 k-i = N 2 k = k. 

Proof. Substituting (2.14) and (2.15) into (2.5) and (2.6) and expand them in series, we 
get the following equations 

A" + cot 9A[ + 2 • 3Ai = -1, (2.17) 

A' 2 ' + cot 04, + 3 • AA 2 = -\ - 6A, - Ub[) 2 + Bf) , (2.18) 

2 sin 6 V / 

B'[ - cot 9B[ + 1-251 = 0, (2.19) 
B 2 - cot 9B' 2 + 2 ■ 3B 2 = 2B ± . (2.20) 

Solve the equations, we find that the regular solutions take the forms (2.16). We check the 
succeeding equations by mathematical induction. 

Assume the forms (2.16) hold for n = 2N — 1 and n = 2N. Since the equations (2.5) 
and (2.6) are homogeneous with respect to r, and the highest derivatives with respect to r 
is linear. Calculation shows that the equations for n = 2N + 1 and n = 2N + 2 take the 
following form 

A'i N+1 + cot9A' 2N+1 + (2N + 2)(2N + 3)A 2N+1 = P 1 - 4(iV+ 1)A 2N , (2.21) 
A 2N+2 + cot 9A' 2N+2 + (2N + 3)(2N + A)A 2N+2 

-4b u [2cot9B 2N+1 + (2N + 1)B 2N+1 ] = P 2 - 2{2N + 3)A 2N+1 , (2.22) 
B'zn+i - cot 9B' 2N+l + (2N + 1)(2JV + 2)B 2N+1 = Q x + 4NB 2N , (2.23) 
B 2N+2 - cot 9B' 2N+2 + {2N + 2)(27V + 3)B 2N+2 = Q 2 + 2{2N + l)B 2N+l , (2.24) 

where (Pk,Qk) are polynomials determined by functions (A 2 j_i, A 2 j, B 2 j-i, B 2 j),(j < N) 
and their derivatives, which all take the forms (2.16). Since cos(2n9) = |(e 2n£ '* + e~ 2n9z ), 
under the constraint of the index of (^) n , we find 

N+l N+l 

Pk = y^ j p kn cos(2n9), Q k = ^ q kn cos(2n9), (2.25) 

n=0 n=0 

with constants ipkmlkn)- Substituting (2.25) into (2.21)-(2.24) and solving the solutions 
with finite and symmetrical conditions, the solutions (^2^+2, B 2N+ i) have the forms (2.16), 
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but (A2N+1, B2N+2) w iH introduce extra terms similar to the second kind Legendre functions, 
which include the following factor 

F(9) = In ( I + ° 0S l ) sin 9 cos 9. (2.26) 
\ 1 — cos J 

F{9) also satisfies the boundary condition and symmetry. Nevertheless, the terms including 
this factor will soon be cancelled at the next step of resolution. Then we prove the regular 
solutions always take the forms (2.16). 

(2.16) means that the Weyl-Lewis-Papapetrou coordinate system (2.1) is an ACMC-00 
one. (2.16) can be transformed into the Legendre polynomials, but the Legendre polynomi- 
als are less convenient to calculate the nonlinear terms than the trigonometric functions. If 
the solutions of (U, W) is determined, V can be solved only by (2.7) with boundary condi- 
tion V(oo,9) = 1. Then the stationary axisymmetrical metric of vacuum is determined in 
principle. We determine the coefficients of the series in the next section. 



III. COEFFICIENTS OF THE SERIES AND EXAMPLES 

Similarly to determine the coefficients of the special functions by recursion, the theorem 
2 implies that, the resolution of (2.5) and (2.6) can be transformed into a problem of solving 
coefficients by recursive relation. For convenience of calculation, we make transformation 
U = exp (|), then (2.5) and (2.6) become 

2e u 

r 2 d 2 .u + 2rd r u + cot 9d g u + d 2 e u + ^-^{drWf + {d e W) 2 \ = 0, (3.1) 

r 2 sin 9 

r 2 d 2 W - cot 9d e W + d 2 e W + (r 2 d r ud r W + d ud B W) = 0. (3.2) 
and (2.14)-(2.16) become 

00 

u 



n=l ^ ' n=l ^ ' 

N n -1 N n 

A n =Y. cos ( 2 ^)> R n = b nk[cos(2k9) - 1], (3.4) 

k=0 k=l 

where N n = p±±] . 

Substituting (3.3) and (3.4) into the original equation (3.1) and (3.2), and expanding 
them in series, we get the relations among the parameters (ak n ,bkn)- (see appendix). In 
each term, the parameters (cifc n , b^n) with the larger indexes always take linear form with 
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bigger coefficients, which correspond to the linear Laplace-like operator in (3.1) and (3.2). 
So it is easy to solve the coefficients. However, the equation are underdetermined for the 
parameters (ofc n , £>&„), and we have the following two sequences of free parameters. Let 

a w = -2, a 3 i = mi, a 52 = m 2 , a 2k +i,k = m k , (3.5) 

11 1 

hi = —^Wu b 3 2 = -^w 2 , b 2 k-i,k = ~2k~ w k' {k = 1,2,3, ■■ ■). (3.6) 

where {mk,Wk} are dimensionless numbers determined by the multipole moments of the 
energy-momentum tensor distribution of the source, aio = —2 is determined by (2.9), and the 
factor is introduced to scale the coefficients in the recursive relation. The free parameters 
(3.5) and (3.6) correspond to 'null data' in papers[23, 24], but the definitions are different. 
The coefficients can be recursively expressed as the polynomials of the free parameters(see 
appendix). Solving the coefficients, we get the series solutions of the metric 

Rs 1 /„ ~. \ ( Re\ 1 n,_ „_ s / RA 

w\ 1 + cos 20 — 
\ r 

5 



u = -2— + -mi(l + 3cos20) ( — ) - -w\(l + cos 20) 
r 3 \ r / 2 



2 4 8 9 / -R \ ^ 

[m 2 cos4fl + (-u? 1 2 + -m 2 )cos2fl + — V + — m 2 ] — + 0(r" 6 ), (3.7) 
7 7 35 do \ r J 



[7=1 + - — + (-mi cos 26 + -m l - - — 

r 2\r ) V 2 6 6 ; V r 

[(-^i 2 - cos26 + 1 - i^ 2 - imr] (^) + 0(r~ 5 ), (3.1 

3 1 , fR„ s 









Ml 






r 







(u> 2 cos 26 + -u> 2 + -iui 

5 5 \ r 



[(-u> 2 + -miwi)cos20+— ^1+2^ miwi + — w 3 ] (-M +0(r- 5 )L (3.9) 
V = 1 + f + (icosM + |) + [(-im, + |) cos 2. - i m , + i] )' + 

Comparing (3.8)-(3.9) with (2.9), we learn the physical meaning of R s and w\, 

R s = 2m, J = -wiR 2 = wiiii 2 . (3-H) 
That is, R s corresponds to the Schwarzschild radius, W\ to the angular momentum. 
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U = exp ( -—) , V = exp( —-](—) sin 2 9 1 , W = 0. 



Now we compare the series solution with some exact metrics, i.e. the Curzon, 
Schwarzschild and Kerr solutions. The Curzon metric is diagonal, 

2 

T J \ T 4 V T 

Expanding (3.12) into Taylor series and then comparing them with (3.8)-(3.9), we find that 
all free parameters vanish 

m k = w k = 0, (V7c). (3.13) 

Evidently, such results is caused by the coordinate system. 

For the Kerr and Schwarzschild solution, under some parameter transformation, in the 
canonical coordinate system (2.4) it becomes 



y. _ 2(R + + R_ + R s cosh u)R s cosh u (3 14) 

~ [(R+ - R-) 2 + R s 2 } cosh 2 uj + 2(R+ + R-)R S coshuj + 4R+R. ' 
[{R+ - i?_) 2 + R 2 S ] cosh 2 c; + 2(R + + R^)R S cosh uj 

+ IrTrZ : 



(3.15) 



w _ [Cg+ - R-? - R 2 S ](R+ + R-+ R s cosh(v) sinhuj 

[(R+ - R-) 2 - R s 2 ] cosh 2 co + AR + R_ ' 1 ' ' 

in which 

R± = V r 2 cosh 2 uj ± 2rm cos 6* cosh uj + m 2 , i? s = 2m, tanhcu = — (3.17) 

m 2 

uj = corresponds to the Schwarzschild solution in the Weyl-Lewis-Papapetrou coordinate 
system. The above formalism is introduced for convenience of the following comparison. 
In contrast the Taylor series of the above functions with (3.8)-(3.9), we get the following 
sequence of the free parameters, 

1 3 2 

mi — 1 OL > ITT''. 



3-11 11-19 



2 



7 


17 2 5-7 4 
— oi — n 

28 2 9 


_ 2? 


17-23 


4 3-7-11 6 


2 14 


« + 2 14 « > 


■ 13 o 


3 ■ 5 3 ■ 19 , 59-67 


a ■ 


Q! H ( 

2 20 2 19 


5 3 


3 4 3-41 , 


— a 6 
32 


w 3 = — —a 77— or 

6 2 12 2 11 



mi = — — — — a — 

914 Ql4 

3 2 ■ 5 ■ 11 ■ 13 

m * = ~ - 2 21 " + " V 9 ~ a ' 2 - ~ 2 20 ^ a4 + ::: 2^ a6 2 21 a§ ' (3 ' 18) 

3 5 o 3 4 3-41 , 3 3 - 7 
wi = -a, w 2 = —a - —a° w 3 = —a + -^~ a 



2 7 32 

11 -13 13-47 , 881 - 3-11-13 7 
Wa = — — — a — — a H — —a ot , 

2 15 2 15 2 15 2 15 

where a = tanhijj < 1. The sequences vanish quite fast. 
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a = corresponds to the Schwarzschild solution. In this case, the 'mass moments' 
rrik 7^ 0, which is caused by the coordinates. So to understand these data we need a unified 
standard. 

Evidently, the above calculation shows that any stationary metric with axisymmetry is 
identical with two sequences {m n ,w n }, in which {m n } is mainly related to the multipole 
moments of mass density, and {w n } to the current distribution. On the contrary, for any 
given sequences {m n , w n } with suitable upper bounds, the series is a solution to the Einstein 
equation of vacuum in the region of its convergence. So the relation between the radius of 
convergence of the series and the values of free parameters {m n ,w n } is important. The 
convergence of the multipole moments serieses were discussed in different context. The 
conditions of convergence for the static solutions were established in [14, 23, 25], and the 
conditions for the stationary ones were given in [24, 26]. Although the results ensure the 
convergence of the above series, how to derive the concrete constraint for free parameters is 
quite complicated. 

In the static case, the radius of convergence may be directly derived from the exact 
solution. The result is enlighten, so we give some analysis. Since W = 0, (3.1) becomes 
Laplace equation, and its general asymptotically flat solution is the Weyl class([2],Ch.l8), 



where P n are the Legendre polynomials of n degree. For any N, the function V can be also 
exactly solved from (2.7) [2]. The Curzon solution is the simplest case with N = 0. In the 

case iV = oo, we have 

Theorem 3. For the series solution (3.19) with N = oo, if the free parameters m n satisfy 
the following condition 



where (C > 0, A > 0, K) are given numbers independent of n, then the series solution (3.19) 
and all its derivatives converge in the region 




n=0 



(3.19) 



m n \ < C(n + l) K \ n , 



(3.20) 



r > r = XR S . 



(3.21) 
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Proof. By the property of the Legendre polynomials, for all n we have 



\Pr,\ < h 



^»(cos. 



d 



-(2n - l)sin0P n _i(cos0) + — P n _ 2 (cos#). 

ad 



(3.22) 
(3.23) 



By (3.23) we get 



^ (C0S ' 



^P„- 2 (cosi 



< (2n-l) + 

< (2n - 1) + (2ra - 5) + (2n - 9) + 



—n(n + 1). 
2 V ; 



(3.24) 



Substituting (3.20), (3.22) and (3.24) into (3.19), we find that the series and its derivatives 
are controlled by, 



u < 



n=l 

oo 



n 



\d rU \ < £5>+i)W 

n=l ^ 

~ oo / 

\deu\ < IB^ 1 )^ 2 

n=l ^ 



n+1 



n+1 



(3.25) 
(3.26) 
(3.27) 



The right hand serieses all converge in the region r > r , so (u, d r u, deu) are absolutely 
convergent in the region r > r , and uniformly absolutely convergent in the region r > r\ 
for any given n > ro- Similarly we can check the results for higher order derivatives. 
Integrating (2.7), we find the radius of convergence of V is also ro. The proof is finished. 

In the stationary case with W ^ 0, the situation becomes more complicated due to the 
nonlinear terms, because (a„fc, 6„fc) are polynomials of (rrij,Wj) of j < However, from 
the recursive relations, we find the following inequalities seems true, 



\a nk \ <CM n (n+l) K , \b nk \ <CM n (n+l) K , V(n, k) 
where (K > 0, C > 0) are constants independent of n, and 

us f I II II I— I \ — ^ 

M n = max \Wj\, \mj\ 2 , \Wj\ 2 j . 



(3.28) 



(3.29) 



In principle, the relation (3.28) can be derived from the results in [24, 26]. If (3.28) holds, 
the series solutions (2.14)-(2.16) will be controlled by 

|Z7|<C £(„ + 1)W^) , \W\<C 1 Y,(n + l) K+l (^) , (3.30) 

n=0 ^ ' "=n ^ ' 



n=0 
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where (C ,Ci) are constants independent of n, and 



A 



lim 




(3.31) 



n 



oo 



Then (3.30) also imply that the series and their derivatives converge in the region r > XR S . 



From (3.13) and (3.18) we find that, in the canonical form of the Weyl-Lewis-Papapetrou 
coordinate system (2.4), the Curzon metric has not 'multipole moment' (where it means the 
free parameters), but the Schwarzschild metric has infinite ones. These results are somewhat 
unnatural and puzzle. Evidently, such results are caused by the coordinate system, although 
the canonical form (2.4) is the most convenient one to solve the metric. So how to extract 
the understandable information from the solutions is also an important problem. 

Similar to the concepts of point charge and dipole in the electromagnetism, an ideal 
explanation for the multipole moments should be expressed in the forms of some conserved 
spatial integrals of the source [21]. However, this ideal is associated with how to define 
the covariant generalized functions for nonlinear differential equations, which may have 
not a general solution for the higher order moments. A realistic explanation for these 
free parameters is to solve them by associating the exterior solution with the interior one, 
the results will endow the parameters with concrete values[20]. However, the influence of 
the coordinate system still exists. In [22], the authors suggested two ways of carrying out 
comparison of approximate and exact solutions: one is calculating the multipole structure of 
the Ernst complex potentials for the solutions, and the second is to generating approximate 
solutions from exact ones by expanding the latter in Taylor series with respect to a small 
parameter. 

To interpret the physical meanings of these free parameters, introducing a standard ruler 
may be a convenient choice. As an approach of first step, we find that the Schwarzschild 
metric is a good ruler, because its properties are the simplest and have been well under- 
stood. By comparing the other solutions with this ruler, we can get some definite and 
understandable meanings of these free parameters. 

The standard exterior Schwarzschild space-time is described by 



IV. INTERPRETATION OF THE MULTIPOLE MOMENTS 



( 




( 



) 



l 



ds 2 



1 



1 



R 



dR 2 - R 2 dB 2 - R 2 sin 2 9 dip 2 . 



(4.1) 
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The transformation between (4.1) and (2.4) reads ["J ] 

r = VR 2 - 2mR + m 2 cos 2 6, cos0= ~ m ) cose u 2 ) 

s/R 2 - 2mR + m 2 cos 2 9 

which is valid in the region i? > ^R s (l + sin 6). In the coordinate system (t, R, 9, <p), the 
line element (2.4) becomes 

ds 2 = U(dt + Wdcp) 2 - V[dR 2 + R(R - R s )d& 2 } - U~^R{R - R s ) sin 2 9 dip 2 , (4.3) 

~ / i? 2 sin 2 9 \ , . 

V = 1 + , DD ^ V: (4.4) 



AR(R - R, 

(t, R, 9, ip) is also an ACMC-oo coordinate system. 

For the Curzon solution (3.12), (U — > 0, V — > oo) corresponds to the surface 

R=-R a (l + ame), ^R S <R<R S , (4.5) 

which is an oblate spheroid. This implies the solution is the analytic extension of the vacuum 
produced by an ellipsoid. However, by V in (4.4), we find the solutions are only valid in the 
region R > R s . Expanding (U, V) in Taylor series, we get the free parameters for Curzon 
solution in the form of Legendre polynomials P n (cos9), 

(iSo P » + i34 P2 -iS P *)(t) 6 + ^ 7) - (46) 

V = 1 + — + — +(P P2) — +(Po R 2 Pa) — + 

R \R J K 12 2J \ R J K 168 28 J \ R J 

41 57 /iC 



(P P 2 Pa) — + 0(r" b ). (4.7) 

V 168 560 ; V R ) 

For the Kerr solution (3.14) and (3.16), the multipole moments in the coordinate system 
(t, R, 9, ip) go as follows 

a2 K-^ Q2 + ^-4li p » + ^(§) 5 + ( r " 7 » + 
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JVp 2 + (1-1 a 2 )P ] (-X + 
24 v 2 24 \ R J 

-a 2 P 2 + (-- 
1 16 V 2 

LV 160 280 ; 
LV 448 560 j 

Comparing the dimensionless coefficients of the term P2(cos0) (jf) 3 in (4.6) and (4.8), 
namely, ^ and \a 2 1 we get a definite concept for the relative deformations of each space- 
time. The bigger the absolute value of the free coefficients, the larger the deformation and 
convection of the gravitating source. In the viewpoint of series, the Kerr solution has not 
any speciality. However, besides (R s = 2m, Wi = \a = -^r), whether the other coefficients 
have some relations with the general physical concepts is unclear. 

The solutions in the metric (4.3) can be directly solved from the following equations 

R(R - R s )d 2 R u + (2R - R s )d R u + cot Qd e u + d 2 e u+ 

2e " [(W + «^, (4,0, 



sin 2 6 LV ' R[R-R S 
R(R — R s )d 2 R W — cot &d e W + d@W + [R(R — R s )d R ud R W + d e ud@W} = 0. (4.11) 

The solution is equivalent to (3.7)-(3.10). 



V. DISCUSSION AND CONCLUSION 



The above procedure provides a simple but effective method to solve the series solution of 
stationary and asymptotically flat metric with axisymmetry to any wanted precision. The 
solution is identical with two sequences of free dimensionless parameters, which correspond to 
the usual multipole moments. The free parameters are determined by the energy-momentum 
distribution of gravitating source. For a wide class of given parameters with suitable upper 
bound, the series solution converges in the region r > R s . For a normal star, we always 
have r ^> R s , so the series provides high precise solutions to the stationary metric with 
axisymmetry. 

To interpret the meanings of free parameters, using dimensionless form and setting up 
a standard ruler are meaningful. Only compared with a simple ruler, we can distinguish 
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the differences between the solutions, and get a clear concept of the free parameters. For a 
given star, to solve the free parameters by associating the exterior metric with the interior 
solution, the results will have concrete physical meanings. However, the matching conditions 
on the surface of the star should be carefully discussed. In [27], we find U G C 1 but V G C°. 
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Appendix 



Substituting (3.3) and (3.4) into the original equation (3.1) and (3.2), and expanding 
them in series, we get the relations among the parameters (a^, bk n ) as follows. 



(R \ 2 ( R x 3 

-^J + (-2o3,i + 6a 30 ) M 

[(6a 4 i + 126 2 x ) cos 29 + (-2a 4 i + 12a 40 ) + 206^ 
[((-8a 6j2 + 14a 5j i) + 4&i,i(3ai,o&i,i + 46 2 ,i)) cos 29 
(-4a 5 , 2 - 2a 5 ,i + 20a 5 , ) + 46 ljl (5a 1 , 6 1 , 1 + 126 2 ,i)] (yj + 0(i 

= |(8fe2i + 2aio6n) (^j + [(166 32 + 206 3 i) + 4a 20 6 n + 4a 10 fe 21 ] (^j + 
[(326 4>2 - 203,161,1 + 1201,063,2) cos 29 + (486 4 , 2 + 3664,1) + 
66i,ia 3 ,o - 803,161,1 + 12ai j0 6 3 ,2 + 802,062,1 + 601,063,1] (^J + 
[((486 5 , 3 + 7265,2) + 1601,064,2 + 4a 3 ,i6 2 ,i + 24a 2 , 6 3 , 2 ) cos 29 + 
(726 5 , 3 + 8865,2 + 56654) + 861,104,0 + 16ai,o6 4 , 2 - 803,163,1 + 24a 2 , 6 3 , 2 + 
12a 2 ,o6 3 ,i + 126 2 ,ia 3 , - 804,161,1 + 8oi, 6 4 ,i] (^j + 0(r^ 6 ) J R s sin 2 9. 

The the solutions of the parameters are given by the following recursive relations 

1 1 2 1 2 

010 = -2, o 20 = 0, a 31 = m 1 , a 30 = -m 1 , 041 = -~u>i, a^o = --w u 

2 2 4 8 2 9 

a 52 = m 2 , a 5 i = -w 1 + -m 2 , a 50 = — w 1 + ^m 2 , 
7 7 35 35 

1 1 2 4 1 2 3 

«62 = --^1^2, CbQl = ~Y0 Wl ~~ 5 WlW2 ' fl60 = ~~ Yq^ 1 ~ I^^ 1 ^ 2 ' ' ' ' 

h 1 h 1 h 1 h 1 1 

6n = -yWl' fo 21 = - jWi, 6 32 = 631 = -w 2 - — w u 

13,131 . , 

642 = -^m x wi - —w 2 , 641 = -ggWi + 7^2 - 24W1W1, (5.3) 

1,11,1 115 

653 = — iu 3 , 652 = u^2 H — w 3 , 651 = miWi wi H iu 2 H w 3 

53 g 3, 52 12 2 9 d ' 51 21 105 15 126 

5 1 1 

36 16 96 

1,3 1 1 15 

6r2 = Wi wvm 2 H mvwo H m\W\ w 2 H if 3 

6 168 56 240 120 36 54 d 

11 1 o 51 1 1 25 1 

6«i = miu> 2 Wi H w 1 m 2 miWi H w 2 H io 3 Wi, ■ ■ ■ 

6 480 420 560 35 45 756 420 



r 



(5.1) 



(5.2) 
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The solutions to (4.10) and (4.11) in the metric (4.3) are given by 
^=l-§ + (iM 1 cos2e + lM 1 )(§) 3 + 

[(jilf! - \wf) cos 20 + Ijtf! - jtifl + 0(r-% (5.4) 

? - 1 + 1 + (I) 2 + - 5* cos2e - S M '> (I)" + 0( - 5 ) + ( 1 " S M ' + 

— V + (-—Mi + — V) cos 46 + (--M 1 + — V) cos 2o] ( ^ J , (5.5) 
64 v 32 64 ; v 8 16 7 J \ R J K J 

W = R s sin 2 6 <| ^ + Wl f^) + 0(r" 5 ) + 



r/ 3 v 3 71 , / i? s \ 3 
too Wi cos 26 H — u> 2 H tfi — + 

LV 16 7 5 80 J v -R y 



27 1 „ 9 , _ 239 5 27 , f 



[( tui H — MiWj + -tu 2 ) cos 26 H H Mi^ H «; 2 1 — >, (5.6) 

LV 64 8 A ' 320 24 20 J V R ' 

where 

Mi = m 1 + -, M 2 = m 2 + ^, M 3 = m 3 + ^^, ■■■ (5.7) 

are somewhat 'pure mass' multipole moments deducted the influence of the coordinates (see 
(3.18)). In the case of Schwarzschild metric, we have (M& = Wk = 0,V7c). 



